Raman superradiance and spin lattice of ultracold atoms in optical cavities 
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We investigate synthesis of a hyperfine spin lattice in an atomic Bose-Einstein condensate, with 
two hyperfine spin components, inside a one-dimensional driven high-finesse optical cavity, using 
off-resonant superradiant Raman scattering. Spatio-temporal evolution of the relative population of 
the hyperfine spin modes is examined by numerically solving the coupled cavity-condensate mean 
field equations in the dispersive regime. We find that beyond a certain threshold, depending on 
the transverse and parallel laser pumps and the atom-atom interactions, Raman superradiance and 
a self-organization of the hyperfine spin components simultaneously occur. Furthermore, the spin 
lattice can be made robust by ramping up the pump intensity in time then turning it off after a 
critical value. 
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PACS numbers: 37.30.-K, 42.65.Dr, 42.50.Pq, 37.10.Vz 

An atomic gas inside a high-finesse optical cavity may 
exhibit self-organization when subject to transverse laser 
pump [H, Q ■ In matter-cavity Quantum Electrodynamics 
(QED) systems, the mechanical effect of the electromag- 
netic fields on the motional states of atoms and phase 
shift effect of atomic motion on the fields induce each 
other mutually in a self-consistent loop. Developments 
in trapping of atomic Bose-Einstein Condensates (BEC) 
in high-finesse optical cavities allow for sufficiently 
strong cavity-condensate coupling to realize nonlinear ef- 
fects, such as bistability, even with cavity photon number 
below unity 0-11] • Quite recently, Dicke Superradiance 
(SR) quantum phase transition (QPT) I9l4llll in a BEC- 
cavity system 



12] has been demonstrated 



Dicke QPT for the single-mode BEC inside a high- 
finesse cavity [IH is characterized by an abrupt increase 
in the number of cavity photons, after a certain thresh- 
old of the pump intensity, which is accompanied by bro- 
ken translational symmetry of the condensate with the 
formation of an optical lattice (lij . Pump-cavity pho- 
ton scattering couples the initial zero-momentum state of 
BEC to a superposition of higher recoil momentum states 
10 . Quite a different scenario happens if a condensate of 
atoms with two different hyperfine sates is pumped by a 



17, \M 



laser field far-detuned from atomic transition 
which case, Raman SR [lil |20| may occur during which 
the hyperfine state of atoms changes. 

There has been much interest on multi-mode atom- 
cavity systems recently, such as bosonic Josephson junc- 
tions inside a single-mode cavity 2l| and spin glas ses of 
single-component BEC in a multi-mode cavit y 1221 . Op- 
tical bistability has been studied in spin-1 (23l. |24| and in 
two-mode BECs d, [25[ . Multi-species systems provide a 
very rich platform for investigations of phase transitions, 
in addition to their practical advantages such as faster 
self-organization with lower threshold, [26j and efficient, 
easily interpretable imaging of correlations in phase tran- 
sition by the cavity field 27|, |28( . 
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FIG. 1. (Color online) Left: Schematic drawing of a BEC in 
one dimensional optical cavity subject to parallel and trans- 
verse laser fields. Cavity has a decay rate of k. Right: The 
doublet of lower levels (b and c) of the BEC atoms are coupled 
by cavity field and laser field via the atomic excited state e 
in Raman scattering scheme. Both laser field and cavity field 
are far detuned from the atomic transition frequencies. 



degenerate hyperfine states, | b), | c) in a one dimensional 
single-mode cavity of frequency uj c as shown in Fig. [T] 
The hyperfine states are coupled to an excited state | e) 
by a transverse pump of frequency loq and the cavity field 
in Raman scheme. The cavity is driven by another laser 
with the same frequency In the dispersive regime 
where the fields are far-detuned from the atomic tran- 
sitions, the Hamiltonian of the system [2(| can be 
written as 
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We consider a condensate of N atoms with two non- 



where i/jj(x, t) (ipj(x, t)) is the annihilation (creation) op- 
erator for a bosonic atom at space-time point (x,t). Vj (x) 
is the external trap potential for the state j — b, c, e, lo oc 
is the frequency of b -f-> c transition and {uij} are the 
interaction strengths of atoms in states i and j. The 
parallel laser field strength is denoted by 771 1 and the an- 
nihilation (creation) operator of the cavity mode is a (a'). 
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For a cavity mode with wave number k, Raman scatter- 
ing Hamiltonian (H^ arnan ) has the form 



H 



Raman 



ih / dx iplho(e- iUot + e iWot )ij} b + H.c. 



-ih 



' dxl pt 90 cos(k X )(a + a^ c + H,.,(2) 



where ho and go are the atom-pump and atom-cavity 
dipole interaction strengths, respectively. Transverse 
pump profile is assumed to be wide enough to take ho 
uniform. Dipole approximation is used for the transverse 
direction. After adiabatically eliminating ip e , under the 
condition of Ao — loo ~ uJbe being larger than the excited 
state linewidth, the Hamiltonian reduces to 



dx ipj'Hipj 



fdx^l^M*) 



V h {x) 



j—b,c " i : j—b : c ' 

where, in a rotating frame defined by the unitary opera- 
tor U = exp — iu>ota)a 1 

rR _ ( 2hh 2 
2m dx 2 \ Ao 
+ (hUo cos 2 {kx)(aa) + a 1 a) + V c {x) + hwbc) + <j~ 
+ hq(a + a') cos(kx)(a~ + a + ) 
— h5 c a/a — ihrju(a — a'). (4) 

Here U = ff 2 /A , rj = h go/^a, S c = w -w c , ct+ = \c)(b\ 
and a~ = \b){c\. 

Early stages of the dynamics is strongly influenced by 
quantum fluctuations that triggers the SR. We consider 
the late time dynamics in which the condensate and field 
variables are assumed to be classical [13] ■ The effect of 
quantum fluctuations are introduced by seeding the cav- 
ity field in numerical simulations. In our case seeding 
is physically performed by the parallel pump that drives 
the cavity. The Heisenberg equations of motion in this 
mean-field regime take the following forms 
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2ir\ cos(kx)a r ip c 



(5) 



ipc = - r 
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- (2hU n cos 2 (kx)\a\ 2 + u cc \ip c \ 2 + u bc \^ b \ 2 ) ^c 
2irjcos(kx)a r tpb (6) 
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dx \ipc\ 2 cos 2 (kx) ) a 
dx cos(kx)(ip*ipb + V'bV'c) + V\\i 



(7) 



where a r is the real part of the cavity field a and n is the 
phenomenological decay rate of the cavity [23]. The Ra- 
man SR is identified by an abrupt increase in the num- 
ber of cavity photons n = \a\ 2 while the QPT in the 
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FIG. 2. (Color online) (top panel) Time dependence of 
strength of effective transverse laser pump -q, dynamics of 
cavity photon number n, (middle panel) total magnetization 
Z, and (bottom panel) spatial and dynamical behavior of 
magnetization density Z(x,t). Cavity photon number and 
magnetization indicate SR at which Z(x, t) undergoes self- 
organization around the time t ~ 1/W- Magnetized and de- 
magnetized domains are not well isolated. Here, in units of 
cj r , 77|| = 1000, ui bc = 1, k = 400, U = -0.5, 5 C = 4800, 
A = -4 x 10 6 and N = 48 x 10 3 . 



atomic system is monitored by the magnetization density 
Z(x,t) — jj(\i(jb(x, t)\ 2 — \tp c (x, t)\ 2 ) which is normalized 
by the total number of atoms N = j dx(\ipb{x,t)\ 2 + 
\ip c (x, t)\ 2 ) so that the total magnetization Z can have 
extremum values 1 and —1 when all atoms are in mode 
b or c, respectively. 

Using second order split step method, and assuming all 
atoms are initially in hyperfine spin state | 6), the mean 
field equations are solved numerically to monitor the dy- 
namics of the system. Fig. [5] shows the time dependence 
of n (top panel) and Z (middle panel) as rj increases 
in time scaled by the recoil frequency to r = hk 2 /2m. 
Pump laser is considered to be red-detuned with Ao < 0, 
and thus Uo,r) < 0. In all plots the absolute value 
of r\ is shown for convenience. Bottom panel of Fig. 
[2] demonstrates the spatio-temporal dynamics of mag- 
netization density, where its translational symmetry is 
broken at the threshold of Raman SR. We consider a 
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condensate of length 2A, A being the wavelength of the 
cavity mode, and hence initially Z{x,t = 0) = 0.5. At 
the SR threshold, atoms change their internal states by 
Raman scattering of the transverse pump into the cav- 
ity field. However not all atoms can equally contribute 
due to spatial modulation of the Raman coupling by the 
cavity mode. Raman couplings in ([5]) and (|6]) change 
both the hyperfine spin and the motional states of the 
atoms. The recoil momentum of the atoms changes by 
±hk. Atoms in in state | c) cause a shift in the cavity res- 
onance due to their spatial overlap with the cavity mode 
by — 2Uq J dx \i(j c \ 2 cos 2 (kx) . The cavity mode builds a 
standing wave potential, 2HUq cos 2 (kx)\a\ 2 , whose min- 
ima traps the atoms in | c) at Xj — jX/2 where j is an 
integer. Overlap of different hyperfine spin states act as 
an atomic polarization grating which is enhanced around 
the sites Xj and stimulates even more Raman scattering 
that completes a self-consistency loop. Accumulation of 
the atoms in | c) around the sites Xj while the atoms 
in | b) around the nodes of the cavity mode, which are 
protected from Raman interactions, results in a spatial 
distillation of the atoms in different hyperfine spin states, 
manifested as a self-organized magnetic lattice with a 
lattice constant A/2. Fig. [5] reveals that magnetized do- 
mains are not well isolated from the de-magnetized ones 
and there is no robust pseudospin lattice structure due 
to too many ringing of the SR pulse. 

We can envision a simple scenario where a robust lat- 
tice of well isolated pseudospin domains can be obtained 
if we just turn the transverse pump off once the magne- 
tization reaches its first minimum. This happens when 
rj = 1.4 uj r . Fig. [3] shows the dynamics of n, Z and 
Z{x,t) for the case where rj increases from to 1.4 u) r 
and then is abruptly turned off. Parallel laser pump re- 
mains on but has no effect. Due to the lack of Raman 
coupling, magnetization Z remains constant after switch- 
ing the transverse field off. 

Finally we briefly discuss the effect of the interac- 
tions and the parallel pump on the threshold of Raman 
SR. Presence of a parallel laser pump plays the role of 
usual seeding of supcrradiant modes in the numerical 
mean field treatment of the SR in the late time classi- 
cal regime. The critical strength of transverse pump r\ c 
for Raman SR to happen is numerically determined when 
Z starts dropping significantly after crossing Z = 0.98. 
Top panels of Fig. [4] indicate that significant increase 
of the parallel laser pump strength would reduce, al- 
beit slightly, T] c , but further increase of m would not 
yield to smaller values of r\ c . In Fig. |4] right panels 
present the same data as in left panels but in logarith- 
mic scale. In the presence of atom-atom interactions, the 
system reaches superradiant phase faster and at a lower 
threshold as the comparison of top and bottom panels 
of Fig. 0] reveals. Bottom panels are without the inter- 
actions. The interaction coefficients Uij = 4irh 2 a,ij/m 
are calculated with the s-wave scattering lengths {a.y} 
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FIG. 3. (Color online) Dynamics of cavity photon number 
n (top panel), magnetization Z (middle panel) and magne- 
tization density Z(x, t) (bottom panel) as the intensity of 
transverse laser field rj is increased from to 1.4 ui r and then 
is switched off. Other parameters are the same as Fig. [2] 
Here, the magnetized domains are well isolated by the de- 
magnetized ones, and a robust magnetic lattice is synthesized. 



taken to be a^c = a c b = 90 a.u. ~ 5.95 x 10 3 A and 
abb = a cc = 106 a.u. « 7.01 x 10~ 3 A, with A = 800 ran. 
This analysis verifies that the parallel pump contributes 
weakly to the threshold of Raman SR. 

We conclude that a BEC with two non-degenerate hy- 
perfine spin components in a high-finesse cavity driven 
by transverse and parallel pumps can exhibit Raman SR 
above a critical value of transverse field strength. Simul- 
taneously, BEC undergoes a phase transition, associated 
both with the external and internal degrees of freedom, 
during which atoms scatter transverse laser field into the 
cavity mode and in return their hyperfine state changes. 
As a result, cavity photon number rises abruptly and at 
the same time there is a sudden increase in the popu- 
lation of higher hyperfine state at periodical positions. 
Self-organization yields spatial distillation of atoms in 
different hyperfine states exhibiting a magnetic lattice 
configuration. 

Even though present analysis is in the mean field 
regime, we can envision that hyperfine spins at different 
lattice sites would be entangled as they interact with the 
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FIG. 4. (Color online) Numerically-obtained critical trans- 
verse laser strength rj c as a function of the strength of parallel 
laser r\\ | , for the cases with (top panels) and without (bottom 
panels) atom-atom interaction. Right panels show same data 
as left panels in logarithmic scale. Atomic scattering lengths 
are given in the text and other parameters are the same as 
Fig- El 



common cavity field, following the resonant entanglement 
of atoms in multitraps scenario [29j . Availability of large 
number of spins per site could make the system advanta- 
geous for explorations of magnetic supersolid properties. 
In contrast to Rayleigh SR, Raman SR can be used for 
a source of entangled photon-spin pairs. Synthesis and 
probing lattice spin models in a single compact cavity 
design, with fast SR induced QPT and self-organization 
properties, promise unique opportunities for quantum in- 
formation applications as well as monitoring QPTs and 
spin correlations with Raman scheme. 
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